Journal of Heredity

Strong population contractions can leave a persistent genomic legacy that can influence populations long after their demographic recovery. While bottlenecks facilitate the removal of strongly deleterious mutations, the effectiveness of purging may be limited in historically small populations. The Kirtlands warbler (Setophaga kirtlandii) is a rare North American songbird with an ancestrally small population. After narrowly evading extinction, they are one of few species that have been delisted from federal protections in the USA. Despite their recovery, a previous study showed evidence for recent inbreeding and a high burden of deleterious mutations that may have not been purged despite strong bottlenecks. Historical DNA offers a unique opportunity to understand the consequences of recent demographic declines on genetic diversity. Here, we use DNA from over 100-year-old museum specimens to estimate changes in genetic load in the Kirtlands warblers pre- and post-bottleneck. We validate our results with forward-in-time genetic simulations and explore how sample size and missing data can affect estimates. Both empirical data and simulations suggest a reduced ability to purge deleterious mutations in this historically small population. Our simulations also highlight that limited sampling design and data quality can constrain the ability to detect changes.

The fitness of immigrants and their descendants determines the effectiveness of gene flow. Genetic incompatibilities or outbreeding depression can limit the spread of novel alleles, while highly fit immigrant lineages can hasten introgression. These fitness effects of gene flow can also differ between generations as immigrant and resident haplotypes recombine. Understanding the genetic factors that shape immigrant fitness over multiple generations is increasingly important as habitat fragmentation threatens populations by reducing genetic variation and leading to increased levels of inbreeding. Few studies have measured the multigenerational fitness of immigrant lineages, especially within populations that had histories of high gene flow. We used 33 years of life history and pedigree data on a population of Florida scrub-jays (Aphelocoma coerulescens) with historically high immigration to quantify the fitness of immigrants and their descendants. We compared the fitness of immigrants and residents as well as their resulting descendants (F1, F2, etc.) to determine the composite genetic effects responsible for fitness differences. We found evidence of additive benefits of immigrant ancestry and heterosis driven by non-additive effects that persists for multiple generations. These results are promising for conservation efforts aiming to increase connectivity and illustrate the complex dynamics that determine the rates of introgression in natural populations.

1Effective population size (Ne) is a critical parameter for evaluating the evolutionary and persistence potential of endangered populations and for designing sustainable conservation strategies. Captive breeding and release programs are widely used across taxa to reduce risk of extinction when natural reproduction is insufficient or no longer possible, making it essential to assess their consequences. We used the case study of the landlocked Saimaa salmon (Salmo salar), one of the most critically en-dangered salmonid populations in Europe, with unique evolutionary significance due to its isolation from other populations since the last glaciation. Using long-term demographic data (1969-2024) from wild-caught founders of a captive breeding and release program, we estimated the effective population size under multiple scenarios of variance in reproductive success. Across scenarios, Ne ranged from 33 to 81 individuals, representing 32%-75% of the census size. Captive breeding practices aimed at equalizing parental contributions during fertilization and early life stages increased Ne by 12% compared to natural reproductive conditions. However, variation in survival after early developmental stages, typically beyond direct management control, remained a key determinant of Ne. Despite recent increases in the number of founders, the population remains genetically vulnerable due to historical bottlenecks. These results highlight that while captive breeding programs can partially mitigate genetic risks, their effectiveness depends critically on both controlled and uncontrolled sources of variance in reproductive success. Strengthening such programs may require combining breeding management with habitat restoration and, where appropriate, genetic rescue to ensure the long-term evolutionary potential of such unique and endangered populations.

Whole genome duplication is a common mutation in eukaryotes with far-reaching phenotypic effects. The resulting morphological, physiological, and fitness consequences and how they affect the survival probability of newly polyploid lineages are intensively studied, but very little is known about the effect of genome doubling on the short-term evolvability of populations. Understanding the effect of polyploidization on the adaptive potential of populations is of crucial importance to predict the future of polyploid populations. In this paper, I investigate the immediate consequences of genome doubling on the genetic variance of populations. To do so, I performed numerical iterations and simulations of how the genetic variance of a quantitative trait changes after polyploidization, under different genetic architectures (additivity, dominance, and epistasis). I found that genetic variance generally decreases after genome doubling. Non-additive gene actions can make autotetraploid populations genetically more diverse than their diploid progenitors in rare cases, notably with overdominance and directional epistasis. By collecting estimates from the agronomic literature, I found that both dominance and epistatic variance contribute to the genetic variance of polyploid populations. These results bring new insights into the adaptive potential of newly formed tetraploid populations, and call for further experimental investigations of how polyploidization is associated with a short-term decrease in evolvability.

The field of genetics of bird migration advances, driven by exponential refinements of sequencing and tracking technologies. In willow warblers (Phylloscopus trochilus), a complex repeat-rich region named MARB (Migration Associated Repeat Block) has recently been found to correlate with the routes taken by individual birds from Europe to their African wintering grounds. However, the genomic location of this region remains unknown. Here, we characterized MARB using a combination of approaches to understand how it evolved. We describe the region using long-read genome assemblies of two willow warbler subspecies (P. t. trochilus and P. t. acredula), two related species, the common chiffchaff (P. collybita) and the greenish warbler (P. trochiloides), and whole genome sequencing data from 76 willow warblers. Finally, we applied karyotyping and fluorescent in situ hybridization techniques on willow warbler spermatocytes to cytogenetically locate MARB. Due to the many repeats, we cannot order scaffolds in silico, but probe hybridization on the karyotype shows that MARB constitutes a single locus (~27.5 Mb) spanning most of the 11th largest chromosome in the willow warbler genome. Interestingly, the MARB regions of all species share several characteristics such as relatively high GC content (50%), a high density of specific repeat families and notably, more than 800 olfactory receptor sequences. Regions homologous to MARB may exist in several migrant bird genomes, though currently unassembled due to their complexity. Resolving these in species with similar migratory polymorphisms to willow warblers will be essential to determine whether MARB influences migratory behaviour across species.

The marbled teal (Marmaronetta angustirostris) is a widely distributed but declining waterfowl species, classified as Near Threatened globally and Critically Endangered in Spain. Despite ongoing conservation actions, including ex situ management and population reinforcement programmes, the genomic consequences of long-term captivity, inbreeding, and patterns of functional genetic variation remain unknown due to the absence of a species-specific reference genome. Here, we present the first chromosome-level genome assembly for this species. The genome was generated using PacBio HiFi long reads and Omni-C data, yielding a 1.15Gb assembly with a scaffold N50 of 76.95Mb. A total of 97.16% of the assembly was anchored into 36 chromosome-scale scaffolds, including the Z and W sex chromosomes. BUSCO analysis recovered 99.2% of conserved avian genes. Gene prediction was performed using both ab initio and homology-based strategies, resulting in 16,048 protein-coding genes. This resource provides a foundation for genomewide analyses of inbreeding, demographic history, and adaptive variation, and will support evidencebased in situ and ex situ conservation strategies for this threatened species.

AO_SCPLOWBSTRACTC_SCPLOWClose-Kin Mark-Recapture (CKMR) is a statistical framework for estimating demographic parameters of wild populations. Instead of recapturing individuals, it relies on the identification of closely-related pairs such as parents and offspring, or siblings. By measuring how often such close-kin are "recaptured" among sampled animals (whether alive or dead), scientists can estimate demographic parameters such as census size, mortality rates, and connectivity. CKMR is starting to change fisheries and wildlife management by giving more reliable demographic information, even for many species that resist conventional approaches. Here we introduce the kinference R package, which provides a set of tools for finding close-kin pairs among thousands of samples each genotyped at thousands of SNPs, and for associated quality control. The CKMR context implies different requirements and assumptions to many other kinship programs. In particular, kinference accounts empirically for linkage without requiring a genome assembly, is able to estimate and control false-negative and false-positive probabilities, and can cope with null alleles. The package has been developed and used in numerous CKMR projects since 2017. This paper documents the assumptions, statistical algorithms, and intended workflow for kinference.

Holcosus orcesi, the Orces Blue Whiptail, is a Critically Endangered lizard endemic to the upper Jubones River basin in southern Ecuador. Restricted to a narrow elevational range within semi-arid Andean shrublands, it represents one of the few montane members of a predominantly lowland lineage. Here we present the first high-quality reference genome for H. orcesi, generated using Oxford Nanopore Technologies long-read sequencing. The assembly spans 1.68 Gb across only 91 contigs, with an N50 of 76.2 Mb and a BUSCO completeness of 96.8%, making it among the most contiguous and complete squamate genomes to date. Structural annotation predicted 25,682 genes, of which 85% showed homology to known proteins and 45% were assigned Gene Ontology terms. Repetitive elements accounted for 46.3% of the genome, with LINEs representing the predominant class. This genome provides a foundational resource for future evolutionary, comparative and conservation-genomic research of H. orcesi and other mountain reptiles, enabling studies of population genomics, local adaptation, and genomic erosion in isolated populations. By expanding the genomic representation of tropical montane reptiles, this work helps address longstanding phylogenetic and geographic gaps in global biodiversity genomics and provides a foundation for evidence-based conservation of H. orcesi and related taxa.

Interest in quantifying linkage disequilibrium (LD, non-random associations of alleles at different loci) has skyrocketed in recent years as researchers have focused on use of LD in genome-wide association studies (GWAS), for studying historical demography, and for estimating effective population size (Ne). The most widely used LD metric is r2 = the squared correlation of alleles at a pair of loci. Despite a half century of efforts, developing an unbiased expectation of r2 as a function of the many factors that can affect it (physical linkage, genetic drift, selection, migration, mutation, mating systems) remains elusive. Furthermore, even when all of these other factors are absent, empirical estimates of r2 are upwardly biased by sampling a finite number (S) of individuals, and that must be accounted for if one wants to focus on the desired signal of LD. Previous approaches to estimate [Formula] have been shown to be biased to greater or lesser degrees. The purpose of this short paper is to demonstrate that a simple and apparently exact expression for [Formula] does exist for the special case where sampling error is the only factor contributing to r2, in which case [Formula] = 1/(S - 1). When other factors contribute heavily to LD, [Formula] shrinks toward 0 as empirical r2 [->] 1. However, for estimating contemporary Ne with unlinked markers, empirical r2 will generally be small and 1/(S - 1) will provide a robust estimate of [Formula].

During meiosis, homologous chromosomes pair to form synaptonemal complexes (SCs) and exchange genetic material through a process known as meiotic recombination. First, programmed DNA double-strand breaks form, followed by the assembly of recombination foci on SCs. These foci mark the sites of recombination intermediates and future crossovers. Distributions of recombination foci along SCs have been studied in many eukaryotes, revealing the interplay between recombination patterns and genome evolution. However, in fish, data on recombination patterns are scarce, and, for the majority of groups, completely absent. Here, we measure the positions of MLH1 foci in 3,504 SCs from 219 male meiotic cells of an African annual killifish Nothobranchius virgatus, a representative of a genus with remarkable karyotype and genome diversity, and present a detailed statistical analysis of its recombination patterns. We found that, in contrast to the several other fish species characterised to date, recombination in N. virgatus occurs across almost entire chromosome arms, excluding (peri)centromeres and telomeres. In the longest SCs, we observed a proximal and a distal peak of the recombination focus frequency and explained the peaks by chromosome pairing dynamics. We also revealed the typical positions of focus pairs, demonstrated interference between foci, with the minimal interfocus distance of 4 m, and described regions of the total recombination suppression near centromeres and telomeres. In sum, our study provides a detailed analysis of recombination patterns in a killifish with a fully acrocentric karyotype and contributes to cytogenomic and statistical methodology for future exploration of meiotic recombination patterns.

O_LIEffective conservation of highly migratory species requires understanding genetic structure across breeding populations and access high{square}resolution markers capable of assigning individuals from mixed aggregates (e.g. bycatch or new nesting sites) to their natal origins. Genomic approaches provide unprecedented resolution but add methodological challenges; thus, it is essential to first build a genomic baseline from known breeding areas and then evaluate strategies for assigning unknown individuals. C_LIO_LITo address this, we used 2b-RAD sequencing, a genomic reduction technique useful for degraded DNA, and loggerhead turtles as a case study. This species shows philopatric breeding, while juveniles and adults form mixed aggregations in foraging grounds. C_LIO_LIOur results highlight the importance of building baselines that include all potential source populations contributing to mixed aggregations. We detected hierarchical genetic differentiation and high resolution and successfully assigned the natal origin of 124 unknown individuals from four Mediterranean foraging grounds. These grounds showed distinct source contributions, and comparisons with previous studies suggest possible temporal shifts in stock composition. C_LIO_LIWe provide a comprehensive genomic baseline for individual assignment of Altanto-Mediterranean loggerhead turtles of unknown natal origin and a general framework for identifying population-specific threats in highly migratory species. C_LI

Social interactions are ubiquitous in nature and have the potential to affect trait evolution, particularly in group-living animals such as cooperative breeders. Interactions among conspecific individuals can affect the amount of additive genetic variation for a trait when the phenotype of an individual is also affected by the genotype of its social partner(s) via indirect genetic effects. Thus, quantifying both direct and indirect genetic effects of social partners is critical for understanding and predicting evolutionary trajectories. While much is known about maternal indirect genetic effects, empirical estimates of indirect genetic effects from other social partners remain limited, particularly in wild populations. Here, we use animal models to assess the contribution of indirect genetic effects from all social partners in a family group (mothers, fathers, and helpers) on juvenile morphometric traits across ontogeny in the cooperatively-breeding Florida scrub-jay (Aphelocoma coerulescens). We found indirect genetic effects of helpers and fathers on nestling weight, but no indirect genetic effect of mothers. Across ontogeny, we found increasing additive genetic variation in both weight and tarsus length. Our study provides a comprehensive assessment of within-group indirect genetic effects in a cooperative breeder and highlights the importance of considering indirect genetic effects beyond maternal effects.

PremiseCentaurea melitensis (Asteraceae) is a problematic invader of grasslands globally, but little is known about its genetic makeup. Here we develop a reference genome to facilitate studies of its invasion history, genetic variation, and evolution. MethodsInbred offspring of a single individual of C. melitensis from its invasion of California, USA were used for flow cytometry to estimate genome size, and for genomic DNA extraction. DNA was sequenced with PacBio HiFi technology (yield = 85.7Gb). The genome was assembled with Hifiasm and annotated with BRAKER3. GENESPACE was used to compare gene order (synteny) with three other species within the subfamily Cichorioideae. ResultsWe estimated a mean genome size of 795.0 Mbp for C. melitensis, and our assembly totaled 696.6 Mbp in 48 contigs (N50 = 55.6 Mbp; BUSCO = 98%), with annotation of 25,157 protein-encoding genes. This included four telomere-to-telomere putative chromosomes, nine additional chromosome arms terminated by telomeric repeats, and a complete chloroplast genome. Synteny varied markedly across the genus and subfamily, suggesting a dynamic history of structural variation in the lineage of C. melitensis. DiscussionWe provide a highly complete and contiguous genome assembly to facilitate the further study of genomic variation in C. melitensis.

AO_SCPLOWBSTRACTC_SCPLOWEnvironmental heterogeneity across freshwater systems often promotes phenotypic variation, yet disentangling environmentally induced variation from heritable differentiation remains a central goal in evolutionary ecology. We investigated the geographic distribution and morphological differentiation, and heritability of shell traits among populations of the freshwater lymnaeid snail Pectinidens diaphanus in Patagonia. Extensive field surveys across 196 freshwater sites revealed that the species occupies a broad range of lentic and lotic habitats and constitutes the only lymnaeid inhabiting southern Patagonia. While reproductive anatomical structures were conserved across populations, shell shape differed markedly among populations from contrasting habitat types, with population identity explaining nearly 50% of total shape variation. Populations from hydrologically unstable habitats (ponds and streams) exhibited more elongated shells and relatively smaller apertures, a pattern consistent with functional responses to hydroperiod variability and desiccation risk. To assess the heritability of this differentiation, we conducted a common-garden experiment across two generations. Shell shape differences between permanent- (lagoon) and temporary- (pond) habitat-derived populations persisted into the G2 generation reared under standardized laboratory conditions, indicating that the observed variation is not solely a response to local environmental conditions but includes a heritable component. Together, our findings demonstrate that P. diaphanus constitutes the sole lymnaeid across southern Patagonia, occupying a broader range than previously documented, and that populations show heritable shell differentiation potentially associated with contrasting freshwater habitats. By integrating large-scale biogeographic surveys with morphometric and experimental approaches, this study provides new insight into how habitat variation may contribute to ecological and evolutionary differentiation in freshwater gastropods.

Prolonged association between mothers and their offspring is common in ungulates, with the level of maternal investment likely to play a central role in shaping this trait. Here we examined patterns of association between mothers and offspring over time, the apparent benefits of association to offspring, and costs to mothers. We analyzed 40 years worth of census data from an individually-monitored, food-limited population of red deer (Cervus elaphus) on the Isle of Rum, Scotland. Starting from birth, female calves associated more frequently with their mothers than male calves in their first year. Calves also associated less with their mothers if the mother did not conceive a new calf. Association frequency decreased with mothers age and population density, and survival over the first year was not related to mother-calf association. Yearlings, now in their second year, were more often associated with their mothers if they were female, if there was no subsequent calf (or the subsequent calf died as a neonate), and if they were still being suckled. Increased association between mothers and yearlings was associated with increased survival to adulthood at 28 months, but suckling a yearling did not improve its probability of survival. For individuals that reached maturity, increased association in the yearling year was associated with slightly shorter adult life spans. The level of association between a calf and mother was not associated with the mothers immediate survival or fecundity. Our findings suggest that juveniles born to poor-condition mothers benefit from prolonged association through improved yearling survival.

In androdioecious species like Caenorhabditis elegans, where the primary mode of reproduction is self-fertilization, the evolutionary role of males has long puzzled biologists. One proposed benefit of males is the potential to escape inbreeding depression. We tested this by enforcing seven generations of inbreeding across nine C. elegans strains differing in baseline male frequency and measuring competitive relative fitness before and after inbreeding. We then relaxed inbreeding for four generations to assess recovery. We predicted that strains with higher male frequency, and greater opportunity for outcrossing, would exhibit faster recovery once inbreeding was relaxed. Strains varied substantially in their responses with most showing significant fitness declines and partial recovery but neither the magnitude of inbreeding depression nor the extent of recovery correlated with male frequency. These results show that male frequency is a poor predictor of inbreeding responses and does not reliably reflect realized outcrossing or its fitness consequences.

Sex ratio distorters (SRDs) are heritable elements that bias offspring sex ratios to enhance their transmission. In the terrestrial isopod Armadillidium vulgare, feminization of genetic males can occur through vertical transmission of the sex ratio distorter known as the f-element, as well as through infection by Wolbachia, a maternally inherited bacterial endosymbiont that can alter host reproduction. Previous studies have focused on the distribution of SRDs and their associations with mitochondrial haplotypes in native European populations, but these patterns are poorly understood in the United States. In this study, we sampled A. vulgare in 12 states, screening individuals for Wolbachia infection, the presence of the f-element, and mitochondrial haplotypes. We found that Wolbachia shows a heterogeneous distribution across populations and haplotypes, in contrast with stronger associations in Europe. The f-element occurred in lower overall frequencies but showed a strong association with mitochondrial haplotype VI. These results indicate that patterns associated with SRD differ from those observed in Europe and suggest that multiple introductions and population mixing have shaped these distributions.

Matching habitat choice provides a mechanism for individuals to maximise their expected fitness by selecting an environment that better fits their phenotype. Many animals choose their local environment by evaluating levels of perceived predation risk against possible resource gain. To test if predation risk is a major driver of habitat choice, we quantify scototaxis, or preference for dark versus light backgrounds, in juvenile guppies. As light backgrounds increase visibility to predators, this aspect of habitat choice captures variation in boldness in small fishes. By rearing and testing 586 fish descended from ten natural populations from Trinidad under common garden conditions, we first quantify (broad sense) heritable variation, i.e. evolutionary potential, within populations. Next, we test for evolutionary divergence among populations in mean preference, and if present, whether ancestral predation regime is a mediator of divergence. Finally, we ask whether families and/or populations differ in the amount of behavioural variation they contain. Habitat choice varied among families (12% of total variance), consistent with heritable variation (0.2). We also found mean preference varies among populations (11% of total variance explained). Evolutionary divergence among-populations is partly explained by ancestral predation regime, with populations from low-predation sites showing a stronger average preference for dark backgrounds than high-predation populations from the same river. Additionally, we find that within-population behavioural variation is greater in high-predation populations. We conclude that guppy populations contain heritable variation that could facilitate adaptive evolution if scototaxis is subject to natural selection. Furthermore, while genetic drift may also contribute to evolutionary divergence among-populations, observed patterns are qualitatively consistent with local adaption to predation regime. Our results suggests that high predation sites favour bolder habitat choice on average, but also that local predation regime shape the evolutionary dynamics of variation, perhaps by maintaining shy-bold variation among-individuals or by favouring individuals with less-predicable behaviour.

In many oviparous animals, egg yolk is the sole source of nutrition until feeding begins, and carbohydrates are present in only small amounts in the yolk. Glucose plays an important role in the developmental processes of various animals. In addition, gluconeogenesis has been reported to occur in the yolk syncytial layer (YSL) of cartilaginous fish and teleosts. In contrast, the role of gluconeogenesis in tetrapods remains unclear. In this study, we used Xenopus tropicalis, an anuran amphibian, which lacks YSL, and therefore provide an opportunity to examine the evolutionary conservation of gluconeogenic mechanisms among vertebrates. In X. tropicalis, liquid chromatography/mass spectrometry revealed that glucose levels increased before liver formation. Subsequent tracer experiments using 13C-labeled metabolic substrates detected gluconeogenesis activity from glycerol and lactate. Expression analyses showed that gluconeogenic genes are expressed in the epidermis and endoderm. Consistently, G0 knockout of fbp1, a key gluconeogenic gene, resulted in a significant reduction in glucose levels, affecting brain development. These findings first demonstrate that gluconeogenesis supports development of X. tropicalis. To the best of our knowledge, gluconeogenesis in developing epidermis has not been reported, highlighting previously unrecognized diversity in tissue-specific metabolism during vertebrate development. Comparative analyses across species will provide further insights into the evolution and functional significance of embryonic gluconeogenesis and nutrient metabolism.

The fungal crop pathogen Blumeria hordei, causal agent of powdery mildew on barley, presents life-history and epidemiological characteristics, as well as and selective pressures due to modern agriculture leading to expected sweepstakes reproduction, that is highly skewed offspring distributions. Using genome-wide polymorphism data and population genomics inferences, we aim to 1) infer the past demographic history and the strength of sweepstakes reproduction in B. hordei, and 2) quantify the contributions of these selective and neutral processes in the genome. An new inference method based on Neural Posterior Estimation and diversity and linkage disequilibrium statistics was developed and tested on simulated and B. hordei genomic data. We confirm that B. hordei exhibits a moderate sweepstakes reproduction (-parameter of 1.6). We highlight that the Site Frequency Spectrum (SFS) appears sensitive to the joint occurrence of sweepstakes and recent demographic changes, which may caution on the reliability of the SFS to infer sweepstakes reproduction. We then scan the genome for selective sweeps, adjusting the significance thresholds of the methods for demographic history and sweepstakes reproduction, thereby yielding a counterintuitive result. When conditioning the significance threshold for sweep detection on simulations under sweepstakes and demography, a very large number of putatively selected regions is found (11.6% of the genome). We suggest that sweepstakes reproduction in B. hordei is due to 1) neutrality (clonal/sexual phases and Boom-and-Bust cycles) generating a genome-wide level of background noise in the coalescent genealogies, and 2) selective sweepstakes due to pervasive positive selection. Our findings have important implications for both population genomic methodology and our understanding of pathogen evolution.